Topographic design in wearable MXene sensors with in-sensor machine learning for full-body avatar reconstruction

Wearable strain sensors that detect joint/muscle strain changes become prevalent at human–machine interfaces for full-body motion monitoring. However, most wearable devices cannot offer customizable opportunities to match the sensor characteristics with specific deformation ranges of joints/muscles, resulting in suboptimal performance. Adequate wearable strain sensor design is highly required to achieve user-designated working windows without sacrificing high sensitivity, accompanied with real-time data processing. Herein, wearable Ti3C2Tx MXene sensor modules are fabricated with in-sensor machine learning (ML) models, either functioning via wireless streaming or edge computing, for full-body motion classifications and avatar reconstruction. Through topographic design on piezoresistive nanolayers, the wearable strain sensor modules exhibited ultrahigh sensitivities within the working windows that meet all joint deformation ranges. By integrating the wearable sensors with a ML chip, an edge sensor module is fabricated, enabling in-sensor reconstruction of high-precision avatar animations that mimic continuous full-body motions with an average avatar determination error of 3.5 cm, without additional computing devices.

always needed to process images/videos in external data terminals 17,18 , proposing the challenges of high bandwidth requirement, hardware expense, and power consumption.
An alternative strategy for high-precision full-body motion monitoring and/or avatar reconstruction involves wearable strain sensors that can mechanically conform to dynamic joint surfaces of the human body, allowing the physiological signals to be collected 1,4,6,8,19,20 . For proprioception purposes, the multi-joint motions of the human body require a set of wearable strain sensors capable of achieving high sensitivities in separate strain ranges 8,21,22 . However, the state-of-the-art strain sensors and existing commercial sensors have limited customizable opportunities to tune the sensors' working windows to match the strain changes of targeted joints/muscles 7,23-25 , leading to erroneous sensing signals and low signal-to-noise ratios. Therefore, to satisfy full-body motion monitoring applications, adequate wearable strain sensor design is highly required to achieve user-designated working windows without sacrificing high sensitivity, accompanied with real-time data processing.
Besides tuning the strain sensor characteristics, an additional challenge is to transmit, store, and process the raw sensor data collected through multiple signal acquisition channels 7,21,26 . One facile approach is through the integration of wireless technology (e.g., Bluetooth) with wearable strain sensors, enabling continuous, realtime streaming of multi-channeled sensor data to an external computing device [26][27][28] . On the other hand, an emerging approach is to process the time-resolved sensor data locally (in-sensor), which can largely reduce communication bandwidths and radio power consumption as well as improve data latency and security 28,29 . To the best of our knowledge, integrating edge computing chip(s) with wearable strain sensors (with customizable sensor characteristics) has not been realized in the literature, especially to monitor and analyze the multijoint and multi-mode movements of the human body.
Herein, wearable sensor modules were fabricated with ML models, either functioning via wireless data streaming or in-sensor edge computing, for full-body motion classifications and personalized avatar reconstruction. Ti 3 C 2 T x MXene nanosheets were specifically adopted for the fabrication of piezoresistive nanolayers due to high electrical conductivity, superior conformability, and ease of processing [30][31][32][33][34] . First, by harnessing the interfacial instability during localized thermal contraction, wrinkle-like topographies were heterogeneously created on the piezoresistive MXene nanolayers, the crack propagation behaviors of which were able to be controlled. Through topographic design, the working windows of resulting strain sensors were managed to be tuned from 6 to 84%, which met the strain ranges of all the joints without sacrificing ultrahigh sensitivities (gauge factor (GF) > 1000). Next, by interfacing wearable strain sensors with Bluetooth chips, the wireless sensor module was able to stream multichanneled strain sensing data continuously, which were input to train an Artificial Neural Networks (ANN) model capable of identifying a variety of full-body motions with 100% classification accuracy. Finally, an edge sensor module was fabricated by integrating wearable strain sensors with an edge computing chip, enabling in-sensor Convolutional Neural Network (CNN) to reconstruct personalized avatar animations with an average avatar determination error of 3.5 cm. In comparison with the wireless sensor module, the edge sensor module avoided wireless data streaming and showed 71% less power consumption for full-body avatar reconstruction.
Polymers are often involved to stabilize the assembled nanostructures 30,43 . In this work, three building block units, singlewalled carbon nanotubes (SWNTs), Ti 3 C 2 T x MXene nanosheets, and polyvinyl alcohol (PVA), were selected for the fabrication of piezoresistive nanolayers. Ti 3 C 2 T x MXene nanosheets were specifically selected over other 2D materials because of their high electrical conductivities, superior mechanical properties, intrinsic hydrophilicity, and ease of processing [30][31][32][33][34] . Detailed characterizations of SWNTs and MXene nanosheets are provided in Supplementary Fig. 1, where the diameter and length of SWNTs were 1-2 nm and 5-30 µm, respectively, and the average diameter of as-exfoliated MXene nanosheets was characterized to be ca. 1000 nm.
Next, the dispersions/solution of MXene nanosheets, SWNTs, and PVA were mixed and then underwent vacuum-assisted filtration to deposit a composite nanolayer on a polyvinylidene fluoride (PVDF) membrane. The thickness of as-filtered MXene/SWNT/PVA nanolayer (abbreviated as ps-MXene nanolayer) was fixed at 400 nm for the rest of this study, and the composition of ps-MXene nanolayer was controllable by adjusting the mass ratios of three building block units in the mixture. Supplementary Fig. 2a-c present the top-down and crosssection scanning electron microscope (SEM) images of a planar ps-MXene nanolayer, where SWNTs were highly dispersed and entangled within the MXene multilayer. Further characterizations of ps-MXene nanolayers (including X-ray diffraction (XRD) and Raman spectra) are provided in Supplementary Fig. 2d, e, respectively. The as-filtered ps-MXene nanolayer was then detached from the PVDF membrane in an ethanol bath, and the freestanding ps-MXene nanolayer was transferred onto a thermally responsive polystyrene (PS) shrink film. In this study, the shrink films were specifically customized to contract only in a uniaxial direction above the glass transition temperature (T g ) of PS (ca. 100°C), which were produced by using the roll-to-roll machine shown in Supplementary Fig. 3. After the thermal contraction at 100°C for 120 s, the width of the uniaxial shrink film was reduced 50%, while its length remained constant (see Supplementary Figs. 4 and 5 for detailed characterizations).
By controlling the thermal contraction region(s), the planar ps-MXene nanolayers (named as M p , M indicates the ps-MXene nanolayer, the subscript p refers to the planar feature) were managed to be deformed into different kinds of M n nanolayers (the subscript n refers to the topographic design). Two categories of M n nanolayers were developed in this work: (1) the ps-MXene nanolayers with homogenous topographies (including M p and M w ) and (2) the ps-MXene nanolayers with heterogeneous topographies (including M p-w-p and M w-p-w ).
In the category of homogeneous topographies, a M p nanolayer was shrunk without any constraints, and the resulting nanolayer was fully covered with periodic wrinkles, which was named as M w nanolayer (the subscript w refers to the wrinkle-like feature, Fig. 1a, (1)). Supplementary Figs. 6 and 7 show the digital photos and SEM images of a M w nanolayer, respectively. In the category of heterogenous topographies, the two ends of a M p nanolayer were fixed (Fig. 1a, (2)), so only the middle region showed periodic wrinkles after thermal contraction, which was named as M p-w-p (the subscript w refers to the wrinkle-like feature, p refers to the planar feature, Fig. 1a, (2)). When the middle part of a M p nanolayer was fixed, the two edge regions were shrunk and displayed periodic wrinkles (named as M w-p-w nanolayer, the subscript w refers to the wrinkle-like feature, p refers to the planar feature, Fig. 1a, (3)). Figure 1b presents the SEM images of M w , M p-w-p , and M w-p-w nanolayers, showing that the generation of ps-MXene micro-wrinkles can be locally controlled. Supplementary Fig. 8 further shows the enlarged SEM image of the transition regions of M p-w-p and M w-p-w nanolayers (between planar and textured regions). The wavelength distribution profiles of M p , M w , M p-w-p , and M w-p-w nanolayers are summarized in Fig. 1c, showing that the periodic wrinkles are constrained within the region(s) experiencing localized thermal contraction. The electrical conductivity of a M p nanolayer (at the MXene/ SWNT/PVA ratio of 85/10/5) was confirmed to be 2479 S cm -1 using a four-point probe, and the electrical resistances of all kinds of wrinklelike ps-MXene textures were provided in Supplementary Fig. 9.

Crack propagation behaviors of M n sensors with homogenous and heterogenous topographies
For the fabrication of wearable strain sensors, all the M n nanolayers (including M p , M w , M p-w-p , and M w-p-w ) were first immersed in a dichloromethane (DCM) bath to dissolve the PS substrates, and the M n nanolayers were detached and became freestanding (see Supplementary Fig. 10 and Experimental Section for fabrication details). As shown in the SEM image in Supplementary Fig. 11 and the atomic force microscopy (AFM) profiles in Supplementary Fig. 12, the wrinkle-like textures of a M w nanolayer were slightly relaxed, and the average wavelength increased from 8.1 to 10.2 µm. Next, the freestanding M n nanolayers were transferred onto VHB TM tapes followed by wiring electrical leads to obtain the corresponding M n sensors (as illustrated in Supplementary Figs. 13 and 14).
Finite element analysis (FEA) was conducted to simulate the strain distribution heatmap of each M n sensor under uniaxial strains. With homogenous topographies, the FEA results of M p and M w sensors were compared. As shown in Fig. 2a, i, when the M p sensor was stretched to 120% of its original length, the M p nanolayer experienced large and concentrated strains (>50%) in its central region. In comparison, in Fig. 2a, ii, when the M w sensor was stretched to 120%, the M w nanolayer experienced attenuated localized strains (<20%), and the localized strains were propagated at the valleys of periodic wrinkles. The crack propagation behaviors of M p and M w sensors during continuous strain loading processes were in situ recorded by using a reflection-contrast microscopy. As recorded in Fig. 2b and Supplementary Movie 1, when the M p sensor was stretched to 120%, long and continuous fractures were generated on its piezoresistive nanolayer. On the other hand, when the M w sensor was stretched to 120%, short and zigzag surface cracks were gradually developed (see Fig. 2b and Supplementary Movie 2). As further shown in the SEM images in Supplementary Figs. 15 and 16, the M w nanolayer showed slower crack propagation than the M p nanolayer under stretching, preventing the conductive pathways of nanolayer from being completely cut off.
With heterogenous topographies, the FEA results of M p-w-p and M w-p-w sensors were compared. As shown in Fig. 2c, iii, the M p-w-p sensor showed the strain distribution profile with hybrid features, where the wrinkle-textured regions experienced attenuated strains (<20%) while the localized strain in the planar region quickly propagated above 50%. Similarly, in Fig. 2c, iv, the M w-p-w sensor presented a region-dependent strain distribution profile. As a result, in Fig. 2d, both M p-w-p and M w-p-w sensors showed region-dependent crack propagation behaviors during the strain loading processes: long/continuous fractures emerged in the planar regions, and short/zigzag cracks propagated in the textured regions. By using ImageJ to quantify the FEA results under 120% stretching, the M w-p-w nanolayer showed an average localized strain of 20%, lower than the M p-w-p nanolayer (32%). Therefore, the M w-p-w nanolayer showed slower crack propagation than the M p-w-p nanolayer  (see SEM images in Supplementary Figs. 17 and 18), enabling M w-p-w sensors to sustain higher uniaxial strain (>50%) than the M p-w-p sensor (<30%) before the conductive pathways of nanolayer being completely cut off. From both FEA simulation and experimental results, the crack propagation behaviors of M n sensors can be controlled by creating homogenous and heterogeneous topographies.
Tunable strain sensing characteristics of M n sensors through topographic designs and stretching directions Three characteristics of a strain sensor are generally evaluated, including (1) sensitivity, (2) linear working window, and (3) maximum working strain. The sensitivity of a strain sensor is normally characterized by GF, as defined in Eqs. 1 and 2, where S ε is the relative resistance change at ε strain, ε denotes the applied strain, R 0 and R ε represent the initial resistance and the resistance under ε strain, respectively. On the other hand, the linear working window of a strain sensor is determined by the strain range where its resistance increased linearly with the applied strain. When the strain sensor reaches its maximum resistance, the maximum working strain is abbreviated as ε max . As the M n nanolayers exhibited directional wrinkle-like textures, the stretching direction showed significant effects on the strain sensing characteristics of a M n sensor. Figure 3a, b present the relative resistance change (S ε )-strain (ε) profiles of M n sensors under two stretching directions (parallel to and perpendicular to the wrinkle axes). Under parallel stretching (Fig. 3a), the linear working windows of M p , M w , M p-w-p , and M w-p-w sensors were characterized to be 3-6%, 8-24%, 25-39%, and 35-50%, respectively, with high GF values of 3400, 1160, 1230, and 1470. Also, the ε max of M p , M w , M p-w-p , and M w-p-w sensors were determined to be 6%, 24%, 39%, and 50%, respectively. In comparison, under perpendicular stretching (Fig. 3b), the M w , M w-p-w , M p-w-p sensors demonstrated narrower linear working windows (with an average strain range of 6%) and lower sensitivities (GF~600). In Supplementary Note 1, detailed FEA simulations and in situ SEM studies were conducted on all M n sensors to investigate the effect of stretching directions on their crack propagation behaviors and strain sensing performance. In short, there are two major advantages of selecting parallel stretching over perpendicular stretching, including (1) higher M n sensors' sensitivities and wider linear working windows and (2) more design opportunities via topographic design. Therefore, we adopted the parallel stretching direction for all M n sensors. Detailed comparisons of strain sensor performance between (1) M p and M w sensors and (2) M p-w-p and M w-p-w sensors are supplemented in Supplementary Note 2.
To further examine the influence of the topographic design on the strain sensing performance, several types of M n sensors were fabricated by (1) varying the areal percentages of wrinkle-like region(s) and (2) changing the distribution of wrinkle-like region(s). First, Supplementary Fig. 19 compares the strain sensing curves of various M w-p-w sensors with different areal percentages of wrinklelike regions. With the areal percentages of wrinkle-like regions increasing from 5 to 75%, the resulting M w-p-w nanolayers exhibited more short/zigzag cracks under parallel stretching.  coverage of short/zigzag cracks prevented the conductive pathways from being completely cut off, leading to wider working windows. Therefore, as the areal percentages increased from 5 to 75%, the linear working window of the resulting M w-p-w sensor increased from 12-14% to 45-69%, respectively (Fig. 3c). Second, several other types of M n sensors were fabricated by varying the positions and distribution of wrinkle-like region(s); the total areal percentage of wrinkle-like region(s) was controlled to be the same at 50%. To evaluate the effect of wrinkle texturing distribution on the strain sensing performance, two kinds of M n sensors, including M p-w and M w-p-w-p , were fabricated. Supplementary Fig. 20 shows the FEA results of M p-w , M p-w-p , M w-p-w , and M w-p-w-p , and the average localized strain of each sensor was quantified by ImageJ on the FEA results. Figure 3d compares the strain sensing curves of M p-w and M p-w-p , which have "one" wrinkle-like region. The M p-w-p nanolayer under 120% stretching showed an average localized strain of 32%, while the M p-w nanolayer illustrated a lower average localized strain of 17%. With lower localized strains, the ε max of a M p-w sensor was characterized to be 37%, which was larger than the one of a M p-w-p sensor (ε max = 24%). Figure 3e compares the strain sensing curves of M w-p-w and M w-p-w-p , which have "two" wrinkle-like regions. The M w-p-w nanolayer under 120% stretching showed an average localized strain of 20% (quantified by ImageJ on the FEA result), while the M w-p-w-p nanolayer illustrated a higher average localized strain of 34%. With lower localized strains, the ε max of a M w-p-w sensor was characterized to be 50%, which was higher than the one of a M w-p-w-p sensor (ε max = 32%). From both FEA results and strain sensing performance, setting the wrinkle-like region(s) at the edge position(s) was able to effectively reduce overall localized strains and increase M n sensors' ε max .
The strain sensing stability of all M n sensors was investigated. Figure 3f demonstrates stable relative resistance changes of a M w sensor under repeated uniaxial strains, as the surface cracks were repeatedly generated at the same spots on the reconfigurable microtextures during reversible stretching processes (in situ recorded in Fig. 3g and Supplementary Movie 3). It is worth to note that the flattened plateaus on the signal peaks in Fig. 3f resulted from the default           Fig. 25.
The mechanical properties of all M n sensors were tested. As shown in Supplementary Fig. 26, all M p , M p-w-p , M w , and M w-p-w sensors showed similar Young's moduli of ca. 150 kPa, which were higher than a bare VHB TM tape (106 kPa). In addition, the hysteresis of a M n sensor (U hysteresis ) was quantified by measuring the maximal signal difference between the stretching and releasing processes, as defined in Eq. 3, where S stretching is the relative resistance change signal, (R-R 0 )/R 0 , of a M n sensor during the stretching process, and S releasing is the relative resistance change of a M n sensor during the relaxation process. Based on Supplementary Fig. 27, the hysteresis values of M p , M p-w-p , M w , and M w-p-w sensors were calculated as 18, 27, 28, and 31, respectively, as the hysteresis of a VHB tape increased with the applied strains of 5%, 15%, 25%, and 40% 44,45 . Fabrication reproducibility is highly critical for wearable sensor applications, as the users can avoid re-calibration of as-fabricated M n sensors and use the in-database S ε -ε profiles as the references 23,46,47 . In this study, the fabrication reproducibility was characterized by calculating the signal error of three M n sensor replicates, as defined in Eq. 4, where σ S ε is the standard deviation of S ε under an applied strain (ε). A smaller signal error indicates that M n sensor replicates are with nearly identical strain sensing characteristics, thus showing higher fabrication reproducibility, vice versa. The large signal error resulted from large S ε variations from three M n sensor replicates. According to Supplementary Movie 1, the M p nanolayer exhibited uncontrolled crack propagation behaviors, showing that the large cracks were randomly generated under strains. In comparison, according to Supplementary Movie 2 and 3, the M w nanolayer demonstrated a more controllable crack propagation fashion, showing that the zigzag cracks were constrained along the valleys of periodic wrinkles and emerged repeatedly under strains. Therefore, with one or more wrinkletextured region(s), the M w , M w-p-w , and M p-w-p sensors (three replicates) demonstrate more consistent strain sensing performance and thus lower signal errors <10% (as shown in Fig. 3h). On the other hand, with only planar nanolayers, the M p sensors (three replicates) exhibited the largest signal errors >50%. The effects of nanolayer composition and thickness on the performance of M n sensors were also investigated. By fixing the nanolayer composition (at the MXene/SWNT/PVA ratio of 85/10/5) and thickness (at 400 nm) ( Supplementary Fig. 28a), the ε max of resulting M n sensors was tuned from 6 to 50%, when the nanolayer topography altered from M p to M w-p-w . In comparison, by fixing the nanolayer thickness (at 400 nm) ( Supplementary Fig. 28b), the ε max only increased from 6 to 16% for the M p sensors, from 24 to 47% for the M p-w-p sensors, from 39 to 65% for the M w sensors, and 50% to 84% for the M w-p-w sensors, when the nanolayer composition varied from 85/10/5 to 65/30/5. On the other hand, by fixing the nanolayer composition (at the MXene/SWNT/ PVA ratio of 65/10/5) ( Supplementary Fig. 28c), the ε max decreased from 16 to 14% for the M p sensors, from 47 to 37% for the M p-w-p sensors, from 65 to 48% for the M w sensors, and 84 to 67% for the M w-p-w sensors, when the nanolayer thickness increased from 400 to 800 nm. Supplementary Note 3 provided more discussions of the effects of nanolayer thicknesses on the M n sensors' topographic design and crack propagation behaviors. Supplementary Fig. 29 summarizes the nanolayer composition, thickness, and topography effects on the ε max of M n sensors. Compared with the strategies of adjusting nanolayer composition and thickness, topographic design served as a more effective approach to tuning the M n sensor characteristics. As further summarized in Fig. 3i, the linear working windows of M n sensors (indices I-XII, see fabrication details in Supplementary Table 1) were able to be tuned from 6 to 84%. In comparison with the state-of-the-art strain sensors (summarized in Fig. 3i), our M n sensors showed userdesignated linear working windows without sacrificing their superior strain sensitivities (GF > 1000) 30,33,43,[48][49][50][51][52][53][54][55][56][57][58][59][60][61][62] .

Wireless sensor module for full-body motion classification
With tunable linear working windows, the M n sensors are suitable to match the strain changes of different body joints, and the proprioception information can be collected with high accuracy to classify, track, and reconstruct full-body motions. As shown in Fig. 4a, by measuring the average strain change ranges of seven joints of a volunteer (including back waist, left/right shoulders, left/right elbows, and left/right knees), we designed and fabricated seven M n sensors with joint-matched linear working windows. In detail, the M p sensor with the linear working window of 3-6% was selected to monitor the back waist bending (average strain change~5%). The M p-w-p , M w , and M w-p-w sensors with the windows of 8-24%, 25-39%, and 35-50% were assigned to monitor the movements of shoulders (average strain change~10%), elbows (average strain change~30%), and knees (average strain change~50%), respectively.
The M n sensors in a single type (e.g., all M p sensors) were not able to collect accurate proprioception signals for full-body motion monitoring. When seven M w-p-w sensors were attached to all the joints (Fig. 4b), the S ε signals collected from back waist were too small (<1.0) and hard to be distinguished from noises. On the other hand, when seven M p sensors were attached to all the joints (Fig. 4c-e), the average strain changes of shoulders, elbows, and knees exceeded the linear working window of M p sensor, leading to erroneous S ε signals. Similarly, using single type of M p-w-p or M w sensors was not able to monitor the high-strain movements of left/ right knees (Supplementary Fig. 30). Therefore, it was necessary to use multi-type M n sensors to collect the proprioception signals of full-body motions. As shown in Fig. 4f, the seven M n sensors were able to successfully record the strain sensing signals in a multichanneled fashion for various full-body motions, including (i) left/ right elbow lifting, (ii) left/right shoulder lifting, (iii) squatting, (iv) stooping, (v) walking, and (vi) running.
To facilely transmit the multi-channeled sensor data to a terminal computing device, the seven M n sensors were interfaced with a Bluetooth chip for the fabrication of a wireless sensor module. The wireless sensor module was composed of seven M n sensors with joint-matched linear working windows, an analog-to-digital converter (ADC) with multi-data-acquisition channels, a microcontroller unit (MCU), and a Bluetooth chip. The equivalent circuit is illustrated in Fig. 4g, and the digital photo is shown in Supplementary Fig. 31. Each M n sensor was connected in series with standard resistors, and the resistor value was specifically selected to be 100 kΩ to ensure wireless data transmission with high accuracy (see Supplementary Note 4, Supplementary Fig. 32, and Methods for more details). By applying a constant voltage of 5.0 V (V in ), the voltage outputs (V i ) were measured by the ADC unit and derived in Eq. 5, where V i is the voltage output of i th sensor channel (i = 1-7), V in is the applied voltage (i.e., 5.0 V), R s is the resistance of the standard resistor (i.e., 100 kΩ), R ADC is the input impedance of ADC unit (i.e., 1 MΩ), and R i is the varying resistance of the M n sensor at ith channel. Afterward, the MCU was programmed with the customized codes to collect the multi-channeled voltage outputs in real time (as recorded in Supplementary Fig. 33), which were then sent out by the Bluetooth chip. Our wireless sensor module demonstrated a data transmission rate (S) at ca. 450 bps (bits per second), which was defined in Eq. 6, where A is the data bits of the Bluetooth chip and t is the data transmission time. As we adopted a commercial Bluetooth HC-06, our wireless sensor module maintained its S value above 400 bps in a broad space (>80 meters) and long-term (>100 h) continuous operation. As shown in Fig. 4h, the wirelessly transmitted M n sensor data was input as the training data (Supplementary Table 2 in GitHub) to train a ML model based on ANN. Detailed ML framework is discussed in Supplementary Note 5. Figure 4i shows the t-distributed Stochastic Neighbor Embedding (t-SNE) scatterplot of multi-channeled M n sensor data in response to six full-body motions (see detailed description in Supplementary Note 6), and six distinct clusters were formed after dimension reduction. Without the use of image/video data, the ANN model was able to achieve a accuracy of 100% for full-body motion classification. The accuracy was defined in Eq. 7 and evaluated by using independent testing data (Supplementary Table 3 in GitHub), where Pi is the determined type of ith full-body motion, and T i is the recorded label of ith full-body motion in testing data.

Edge sensor module for in-sensor avatar reconstruction
Although the wireless sensor module ensured facile data connectivity to a terminal computing device, continuous transmission of multichanneled sensor data led to high energy consumption for long-term monitoring 29,63 . In addition, the wireless transmission approach always faces interruption and disconnection challenges, especially when the wireless sensor modules are applied in distant fields or in underwater scenarios 64,65 . To address these challenges, edge data computing in a ML chip becomes an emerging approach to monitoring and determining full-body motions in high accuracy. By integrating seven M n sensors with a ML chip, an edge sensor module was fabricated with the equivalent circuit illustrated in Fig. 5a; the digital photo is shown in Supplementary Fig. 34. The edge sensor module was composed of seven M n sensors with joint-matched linear working windows, a multichanneled ADC unit, a ML chip (with integrated Bluetooth function). It should be noted that the Bluetooth function was used to intermittently send the processed results out of the ML chip. Our edge sensor module with in-sensor ML models was able to reconstruct personalized avatars that mimicked the continuous fullbody motions in high precision and accuracy. Precise monitoring of continuous human motions has long been a primary center for various applications, such as gesture/gait recognition 26,27 , animation production 66,67 , remote healthcare 6,68 , and virtual reality 9,11 . During the coronavirus pandemic, the technologies for sensing delicate body motions (e.g., trembling, shivering) become necessary, as the doctors can monitor the patients' symptoms in real time 1,69 . Also, the human motion sensing technologies have significant impacts in various industries, including sports 4,5 , healthcare 69 , and gaming entertainment 70,71 .
In this work, the in-sensor ML model for avatar reconstruction was built up in two steps. First, the figure aspects and 15 joint locations of a volunteer were extracted from a pre-recorded video through an OPEN POSE program 72 . The detailed description of OPEN POSE program is provided in Supplementary Note 7. As shown in Fig. 5b, the stationary stickman avatar with 15 joints was constructed. Afterward, seven M n sensors with joint-match working windows were attached onto the joints of a volunteer, and the time-resolved M n sensor data were collected in a multi-channeled fashion, which were input as the training data (Supplementary Table 4 in GitHub) to train a ML model offline based on CNN (see details in Supplementary Note 8). The CNN model with optimal hyperparameters was then uploaded to activate the edge sensor module. When the volunteer performed a series of full-body motions, the seven M n sensors were able to monitor the localized strain changes of different joints. As shown in Fig. 5c, the multichanneled resistance change profiles were collected and transformed into the voltage outputs through ADC. By continuously receiving the voltage data, the ML chip with in-sensor CNN model enabled real-time and high-accuracy determination of 15 avatar joint locations, which were then transformed into personalized avatar animations by the OPEN POSE program (as shown in Supplementary Movie 4).
The determination error of in-sensor avatar reconstruction was calculated by measuring the difference between the real joint locations (extracted from camera-recorded videos, Supplementary Table 5 in GitHub) and the determined joint locations (computed by edge sensor module, Supplementary Table 6 in GitHub). The average determination error is defined in Eq. 8, where P i t is the CNN-determined location of i th joint at the time t, T i t is the camera-recorded locations of i th joint at the time t, 170 cm is the volunteer's physical height, and 830 is the corresponding avatar height in the virtual coordinate system, and 15 is the number of monitored joints. By comparing 15 joint locations ( Fig. 5d and Supplementary  Fig. 35), the average avatar determination error was calculated to be 3.5 cm. As shown in Fig. 5e, the CNN-determined avatar animations successfully mimicked the volunteer's full-body motions in high precision and accuracy. It is worth mentioning that the avatar animation sometimes moved ahead of the full-body motions, specifically the squatting motions. As shown in Supplementary Fig. 36 and Supplementary Movie 4, the avatar's squatting movement (at the 22.6th second) was ahead of the video-recorded squatting motion (at the 23.0th second). The ahead motion determination was due to the early signals from the M p sensor on the back waist, where the preparatory muscular stretching happened before the squatting motions (see details in Supplementary Note 9). Supplementary Note 10 and 11 discuss more potentials of the edge sensor modules, including the improvement of the comfort levels of wearable M n sensor modules and the method to re-construct a 3D avatar ( Fig. 5f and Supplementary Movie 5).
The edge sensor module demonstrated a critical advantage in power consumption over the wireless sensor module. Figure 5g compares the power consumption of wireless and edge sensor modules to reconstruct avatar animations (see detailed calculation in Methods). The wireless sensor module showed power consumption of 31.5 mW, the majority of which was used to wirelessly transmit multi-channeled sensor data to a terminal computing device (20 mW). It is worth noting that the power consumption of avatar reconstruction in the terminal computing device was neglected. On the other hand, the edge sensor module consumed 9 mW to achieve the same task of avatar reconstruction, which was 71% lower than the wireless one. The edge sensor module with high power efficiency is beneficial for full-body motion monitoring, as miniature batteries with limited energy capacities are normally used for wearable applications.

Discussion
In this work, through harnessing the interfacial instability during localized thermal contraction, a variety of M n sensors with engineered MXene microtextures were fabricated, where the in-plane crack propagations and thus the strain sensing characteristics were systematically tuned (Fig. 6a). Four kinds of M n sensors with ultrahigh sensitivity (GF > 1000) and joint-matched working windows were attached onto the multi-joint surfaces of a volunteer to monitor fullbody motions with high signal-to-noise ratios. As shown in Fig. 6b, by interfacing M n sensors with a Bluetooth chip, a wireless sensor module was assembled to achieve continuous and real-time streaming of multi-   channeled M n sensor data, and an ANN model was trained with 100% accuracy to classify various full-body motions. By further integrating M n sensors with a ML chip, an edge sensor module was fabricated for the in-sensor reconstruction of personalized avatar animations that mimic diverse full-body motions with an average avatar determination error of 3.5 cm, without external computing devices. Also, the standalone edge sensor module avoided wireless data streaming and demonstrated 71% less power consumption for avatar reconstruction than the wireless sensor module. Our work has demonstrated multiple interdisciplinary advances. Our first advance is to develop a low-cost, scalable, and controllable approach to engineering wrinkle-like MXene microstructures via localized thermal contraction, which cannot be easily achieved by conventional bucking methods (Supplementary Figs. 37 and 38). Our second achievement is to control the crack propagation behaviors of piezoresistive nanolayers through topographic design including the M n nanolayers under different stretching directions, the M n nanolayers with varying the areal percentages of wrinkle-like region(s), and the M n nanolayers with varying positions and distribution of wrinkle-like region(s). As a result, the strain sensing characteristics of wearable strain sensor can be customized. Most of the reported works in the literature focuses on adjusting the piezoresistive nanolayer compositions to pursue high GFs, yet few strategies were developed to customize the linear working windows 23,30,35,47,62 . In this work, the linear working windows of our M n sensors can be tuned without sacrificing their ultrahigh strain sensitivity (GF > 1000). Third, compared with the state-of-the-art works in Table 1, our edge sensor module showed advances in terms of versatile integration method. As mentioned by multiple important review articles 29,69,73,74 , an urgent challenge in wearable sensors is to enable efficient transmission of large amounts of sensor data followed by in-sensor machine learning with power-efficient local computing. Multidisciplinary efforts are required to address this challenge by advancing hardware/software development and optimizing the sensor/ circuit interfaces. In this work, our edge sensor module was carefully designed to achieve full-body motion monitoring coupled with edge data processing and in-sensor machine learning, which has not been reported before. We believe that the integration method (i.e., wearable/stretchable sensors with edge computing chip(s)) demonstrated in this work is highly versatile and can benefit the advances in other fields, including wearable performance devices in sports and underwater soft robots. purification system (Merck Millipore) and used as the water source throughout the work.

Preparation of Ti 3 C 2 T x MXene nanosheets
Ti 3 C 2 T x MXene nanosheets were prepared according to the literature 75 .
1.0 g of LiF was added to 6.0 M HCl solution (20 mL) under vigorous stirring. After the dissolution of LiF, 1.0 g of Ti 3 AlC 2 MAX powder was slowly added into the HF-containing solution. The mixture was kept at 35°C for 24 h. Afterward, the solid residue was washed with DI water several times until the pH value increased to ca. 7.0. Subsequently, the washed residue was added into 100 mL of DI water, ultrasonicated for 1 h, and centrifuged at 1308 g for 30 min. The supernatant was collected as the final suspension of Ti 3 C 2 T x MXene nanosheets with the concentration of ca. 5 mg mL -1 .

Preparation of SWNT dispersion
The SWNT dispersion was obtained by adding the SWNT powders into the SDS solution (at the SDS concentration of 2 mg mL -1 ) at the mass ratio of SWNT:SDS = 1:20. Then, the mixture was ultrasonicated for 2 h by a probe sonicator, and the concentration of the final SWNT dispersion was about 0.1 mg mL -1 .

Calculation of dimension contraction of uniaxial shrink films
A uniaxial shrink film was cut into multiple rectangle-shaped pieces followed by thermal contraction in an oven at 100°C. The dimensions of the uniaxial shrink film before and after thermal contraction were quantified by taking photographs and sampling their gray-scale line profiles using ImageJ. The dimension contraction of a uniaxial shrink film was calculated by Eq. 9, Dimension contraction À ð Þ= 1 À Dimension after thermal contraction Dimension before thermal contraction ð9Þ

Fabrication of freestanding M p nanolayers
Two dispersions of SWNTs and MXene nanosheets were mixed at different ratios, and 5 wt.% of PVA was then added into the SWNT-MXene dispersions. The as-prepared SWNT-MXene-PVA dispersions were next deposited onto PVDF membranes (0.22 µm pore, Merck Millipore) through vacuum-assisted filtration systems. To remove SDS residues, the filtered SWNT-MXene-PVA thin films (abbreviated as ps-MXene nanolayers) were rinsed with excessive DI water. Afterward, freestanding ps-MXene nanolayers were detached from the PVDF membranes by immersing them in ethanol. The planar ps-MXene nanolayers were categorized and abbreviated as M p nanolayers.

Fabrication of freestanding M w nanolayers
A uniaxial shrink film was cut into multiple rectangle-shaped pieces (4 × 8 cm 2 ), washed with ethanol, and dried under N 2 flow. The cut shrink films were next treated with oxygen plasma for 2 min to enhance the hydrophilic interactions between PS substrates and M p nanolayers 76 . Afterward, the M p nanolayers were carefully transferred onto the plasma-treated uniaxial shrink films followed by overnight drying. The M p nanolayer-coated shrink films were then heated in an oven at 100°C for 120 s to induce uniaxial thermal contraction. By harnessing surface instability during thermal contraction, the M p nanolayers were deformed into wrinkle-like microtextures (abbreviated as M w nanolayers). The shrunk samples were then immersed in DCM to dissolve the PS substrates to obtain freestanding M w nanolayers (see SEM image in Fig. 1b), which were sequentially rinsed with DCM, acetone, and ethanol. The M w nanolayers were stored in ethanol.      Fabrication of freestanding M p-w-p nanolayers A uniaxial shrink film was cut into multiple rectangle-shaped pieces (4 × 8 cm 2 ), washed with ethanol, and dried under N 2 flow. The cut shrink films were next treated with oxygen plasma for 2 min to enhance the hydrophilic interactions between PS substrates and M p nanolayers. Afterward, the M p nanolayers were carefully transferred onto the plasma-treated shrink films followed by overnight drying. Then, the two ends of the uniaxial shrink films were fixed by adhering thin glasses at their backside (see schematic illustration in Fig. 1a). The M p nanolayer-coated shrink films were then heated in an oven at 100°C for 120 s to induce uniaxial thermal contraction. By harnessing surface instability during thermal contraction, the M p nanolayers were deformed into the M n nanolayers with wrinkle-like microtextures localized only in the middle region (abbreviated as M p-w-p nanolayer). The shrunk samples were then immersed in DCM to dissolve the PS substrates to obtain freestanding M p-w-p nanolayers (see SEM image in Fig. 1b). It is worth to note that the M p-w-p nanolayers were attached on PS shrink films, and the glass slides were adhered at the backside of shrink film. In other words, the textured M p-w-p nanolayer and the glass slides were attached to the two sides of a shrink film (see Supplementary Fig. 10). After the M p-w-p -coated PS samples were immersed in a DCM bath, the middle PS shrink film was dissolved, and the M p-w-p nanolayers were detached and became freestanding. The freestanding M p-w-p nanolayers were sequentially rinsed with DCM, acetone, and ethanol. The M p-w-p nanolayers were stored in ethanol.

Fabrication of freestanding M w-p-w nanolayers
A uniaxial shrink film was cut into multiple rectangle-shaped pieces (4 × 8 cm 2 ), washed with ethanol, and dried under N 2 flow. The cut shrink films were next treated with oxygen plasma for 2 min to enhance the hydrophilic interactions between PS substrates and M p nanolayers. Afterward, the M p nanolayers were carefully transferred onto the plasma-treated shrink films followed by overnight drying. Then, the middle sections of the uniaxial shrink films were fixed by adhering thin glasses at their backside (see schematic illustration in Fig. 1a). The ps-MXene-coated PS device was then heated in an oven at 100°C for 120 s to induce uniaxial thermal contraction. By harnessing surface instability during thermal contraction, the M p nanolayer was deformed into the M n nanolayers with wrinkle-like microtextures localized only in the edge regions (abbreviated as M w-p-w nanolayer). The shrunk samples were then immersed in DCM to dissolve the PS substrates to obtain freestanding M p-w-p nanolayers (see SEM image in Fig. 1b). It is worth to note that the M w-p-w nanolayers were attached on PS shrink films, and the glass slides were adhered at the backside of shrink film. In other words, the textured M w-p-w nanolayer and the glass slides were attached to the two sides of a shrink film (see Supplementary Fig. 10). After the M w-p-w -coated PS samples were immersed in a DCM bath, the middle PS shrink film was dissolved, and the M p-w-p nanolayers were detached and became freestanding. The freestanding M w-p-w nanolayers were sequentially rinsed with DCM, acetone, and ethanol. The M w-p-w nanolayers were stored in ethanol.

Finite element analysis (FEA) simulation
The 3D models of M p , M w , M p-w-p , and M w-p-w microstructures were built by using SolidWorks 2018, and their surface characteristics were modelled by the Freeform feature. The FEA models of M p , M w , M p-w-p , and M w-p-w nanolayers were further built by the Static Structural module of ANSYS Workbench 19.0 (see Supplementary Fig. 39). The simulation parameters for ps-MXene nanolayers were set as follows: Young's modulus~1.7 GPa, Poisson's ratio 0.227, and mass density 1.25 g cm -3 . Cartesian coordinate was chosen for the mesh method, and the element size was set to be 100 µm. The FEA simulation was conducted as shown in Supplementary Fig. 40, where the left boundary of the 3D model was fixed, while the right boundary was set to be movable along x and y directions. Uniaxial stretching was simulated by moving the right boundary, and the equivalent elastic strains and the overall deformation were recorded.
Fabrication of M p , M w , M p-w-p , and M w-p-w strain sensors (M n sensors) The freestanding M p , M w , M p-w-p , and M w-p-w nanolayers in ethanol were carefully transferred onto VHB TM tapes followed by overnight drying. Copper wires were connected to the two ends of the M n nanolayers, and silver paste was applied between ps-MXene nanolayers and copper wires to ensure good electrical contacts. The resistance profiles of M p , M w , M p-w-p , and M w-p-w strain sensors were monitored by Industrial Multimeters (EX503).
Calculation of crack-to-width ratios (φ) and crack densities (ρ) of M n sensors The crack-to-width ratios (φ) and crack densities (ρ) of M n sensor were investigated via in situ SEM. The definitions of φ and ρ are described in Eqs. 10 and 11, respectively, where max L crack À Á is the length of the longest surface crack, P ðL crack Þ is the cumulative length of all surface cracks, and A and W are the area and width of a M n nanolayer, respectively.

Circuit design of a wireless sensor module
The equivalent circuit of the wireless sensor module is described in Fig. 4g. The wireless sensor module was designed with seven data collection channels, which consisted of seven M n sensors (including one M p , two M p-w-p , two M w , and two M w-p-w sensors), an ADC (AD7606, Risym), a MCU (MCU-PCA9658, DeXin Electronics), and a Bluetooth module (HC-06, XinTaiWei Electronics). Each M n sensor was connected to a standard resistor (100 kΩ) in series. By applying a constant voltage of 5.0 V, the ADC unit measured the voltage outputs (V i s ) in real time, which were derived in Eq. 5. Afterward, the multi-channeled voltage signals were processed by MCU, and sent out by Bluetooth,

Circuit optimization of a wireless sensor module
First, to achieve high sensitivity of wireless sensor module, a highresolution 16-bit ADC (AD7606, Risym) was adopted. Second, standard resistors were used as the regulator to enable high measurement accuracy with low wireless transmission errors, which were calculated by following Eq. 12,

W ireless Transmission Error
where R input is the resistance value of a M n sensor measured directly by a digital multimeter, and R output is the output resistance value of a M n sensor wirelessly transmitted and converted from the ADC-recorded voltage value.
In this work, we measured wireless transmission error by using the M n sensors under different strains. As shown in Supplementary Fig. 32, with 100 kΩ standard resistors, the wireless sensor module exhibited low and strain-stable transmission errors <4%. In the other words, a M wp-w sensor under 48% strain (with a 39.1 kΩ resistance) was read as 37.7 kΩ after wireless data transmission, showing a transmission error of 3.6%. On the other hand, with 5 kΩ standard resistors, the wireless sensor module exhibited higher and strain-dependent transmission errors (the average error >20%). For instance, a M w-p-w sensor under 48% strain (with a 39.1 kΩ resistance) was read as 13.3 kΩ after wireless data transmission, showing a transmission error of 65%.
On the other hand, to simultaneously record the signals of multiple M n sensors, the execute codes were programmed into MCU (the codes were provided in GitHub: https://github.com/jiali1025/ Wearable-MXene-Sensors-with-In-Sensor-Machine-Learning-for-Full-Body-Avatar-Reconstruction).

Power consumption of a wireless sensor module
The power consumption (P) of the wireless sensor module is calculated based on Eqs. 13 and 14 where R s is the resistance of the standard resistor (i.e., 100 kΩ), and R M p , R M pÀwÀp , R M w , and R M wÀpÀw are the resistances of M p , M p-w-p , M w , and M w-p-w sensors, respectively; R module is the resistance of the wireless sensor module, P ADC , P MCU , and P BLE are the power consumptions of ADC (ca. 5 mW), MCU (ca. 5 mW), and Bluetooth units (ca. 20 mW), respectively. V is the applied voltage (i.e., 5.0 V).

Circuit design of an edge sensor module
As shown in Supplementary Fig. 34, seven M n sensors were integrated with an ADC unit and a ML chip (ARDUINO, Element14 Pte Ltd) in series. The equivalent circuit is illustrated in Fig. 5a. Basically, the ADC unit collected multi-channeled M n sensor data and send them to the ML chip. Then, the ML model (which was trained offline) was able to process the M n sensor data followed by the determination of full-body avatar joint locations. The generated avatar joint locations were then transmitted by using the Bluetooth chip integrated in the ML chip. The execute codes were programmed into the ML chip (the codes were provided in GitHub: https://github.com/jiali1025/Wearable-MXene-Sensors-with-In-Sensor-Machine-Learning-for-Full-Body-Avatar-Reconstruction).

Power consumption of an edge sensor module
The power consumption (P) of the edge sensor module is calculated based on Eqs. 15 and 16 P = V R module + P ADC + P MLChip ð15Þ where R module is the resistance of the edge sensor module calculated by Eq. 13, P ADC is the power values of ADC (ca. 5 mW), and V is the applied voltage (i.e., 5 V), I Chip is the current required for the ML chip (the average current dissipation is 21.5 mA). For the ML chip, it required 1,477 regressions for 1 min avatar motion reconstruction, and each regression took about 1 ms.
Characterization XRD analyses were conducted using an X-ray diffractometer (Bruker, D8 Advance X-ray Powder Diffractometer, Cu Kα (λ = 0.154 radiation)) at a scan rate of 4°min -1 . The as-prepared MXene nanosheets were characterized by using a high-resolution transmission electron microscopy (HRTEM, JEOL 2010F). The surface morphologies of M p , M w , M p-w-p , and M w-p-w nanolayers were characterized by using a SEM (FEI Quanta 600) and a field emission SEM (JEOL-JSM-6610LV) operating at 15.0 kV. Surface roughness of M p , M w , M p-w-p , and M w-p-w nanolayers was measured by AFM (Model: Bruker Dimension ICON). The characteristic crack lengths and the crack densities of M p , M w , M pw-p , and M w-p-w nanolayers were quantified by using ImageJ. Fatigue tests were performed on the M w strain sensor for 2000 cycles under repeated uniaxial strains, which were performed on a tensile tester (Instron 5543, Instron, USA) with a 500 N load cell.

Data availability
The data generated in this study are provided in the Supplementary Information/Source Data file. Supporting files of Supplementary Tables 2-6 generated in this study have been deposited in the public GitHub (https://github.com/Haitao008/Supporting-Tables) and Zenodo 77 (https://doi.org/10.5281/zenodo.7012006) without any restrictions. Source data are provided with this paper.